Large thermal Hall coefficient in bismuth 
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We present a systematical study of thermal Hall effect on a bismuth single crystal by measuring 
resistivity, Hall coefficient, and thermal conductivity under magnetic field, which shows a large 
thermal Hall coefficient comparable to the largest one in a semiconductor HgSe. We discuss that this 
is mainly due to a large mobility and a low thermal conductivity comparing theoretical calculations, 
which will give a route for controlling heat current in electronic devices. 



Thermal Hall effect (Righi-Leduc effect) is the thermal 
analog of the Hall effect. A temperature gradient -^ in- 
troduced along x direction yields a temperature gradient 
—r 11 along y direction in the presence of the magnetic 
field B z along z direction. A ratio of the temperature 
gradients is denned as a thermal Hall coefficient -Rth ; 



-Rth = ~^"^f~ — vRh = Mh, 



(1) 



which shows a measure of oriented heat current per mag- 
netic field. The definition of -Rth differs from that of 
the Hall coefficient i?n- ^?th and Rh are inseparably 
related, since they are both governed by the Fermi sur- 
face topology. The relationship between Rth and -Rh (a: 
electrical conductivity) shown in Eq. (1) was derived by 
Bridgman [l| by considering heat and electric transport 
in terms of reversible processes in thermodynamics, and 
qualitatively explains the thermal Hall effect in conven- 
tional metals [2, 3] . The thermal Hall effect of semicon- 
ductors has been studied until 1960s 0, [5J. As a result, 
it is known that -Rth of the semiconductors satisfies a 
relationship |Q|, 



Rth 



^tot ' 



-m 



(2) 



that takes into account heat conduction along the lon- 
gitudinal direction through lattices besides electrons [7|, 
where K e \, K to t(= ^ph + Keh K ph^ lattice thermal con- 
ductivity) , and /ijj represent electron thermal conductiv- 
ity, total thermal conductivity, and Hall mobility, respec- 
tively. Recent developments of this research field have re- 
vealed phonon and magnon thermal Hall effects, in which 
carriers of heat currents along the y direction are not 
electrons but phonons and magnons [gj. These extreme 
effects have a potential for controlling heat current in 
electronic devices as well as a thermal rectifier [9h12i|. 

Bismuth is " an old material" , which has been exten- 
sively studied owing to its unique physical properties such 



'Electronic address: kobayashi.wataru.gfflu. tsukuba. ac.jp 



as large positive magnetoresistance, huge orbital dia- 
magnetism, and peculiar magnetic-field dependent ther- 
moelectric properties [13j |. Recent increasing interests 
in Dirac fermions and strong spin-orbit interaction in a 
heavy element such as Bi, Pt, or Au motivate many re- 
searchers to study this material. Since the bismuth ex- 
hibits large /zh, relatively large k i, and low lattice ther- 
mal conductivity n p h due to both a semimctallic state 
and the heavy mass of Bi element, a large Rth is ex- 
pected through Eqs. (1) and (2). Though Rth of bis- 
muth has been already reported in 1910s [14J, there exists 
the data only at around room temperature and low mag- 
netic field. Thus, in this letter, we report thermal Hall 
effect on a bismuth single crystal up to 3 T and in the 
temperature range of 75 to 300 K, and first analyze tem- 
perature dependent Rth using these equations. We also 
discuss a possible way to further increase -Rth- 

Resistivity, Hall coefficient and thermal conductiv- 
ity measurements of a high quality single crystal (Bi 
99.999%) with a dimension of 1.15 x 2.35 x 5.15 mm 3 were 
performed. For the measurements, electric current i x and 
heat current j x were applied along the bisectrix axis un- 
der the magnetic field introduced parallel to the trigonal 
axis of the bismuth single crystal. Induced Hall voltage 

AT 

Vy and thermal Hall temperature gradient -^ were de- 
tected along the binary axis. The resistivity and Hall 
coefficient were measured by a conventional four probe 
method. 

The thermal conductivity and the thermal Hall coef- 
ficient were measured by a steady state technique. One 
edge of the crystal was glued to a cold head of a closed 
cycle refrigerator by varnish (GE7031), while a resistive 
heater (~ 120 fl, KYOWA strain gauges) was attached on 
the opposite edge. Using the resistive heater, a longitu- 
dinal temperature gradient ^- was generated. Through 
the measurement, the heat delivered to the hot end of 
the crystal was kept constant (1.32W). In the presence of 
the magnetic field (CRYOGENIC, MINI-CFM) applied 
perpendicular to ^p-, the transversal temperature gra- 
dient -£*- was detected along the binary axis. (See Fig. 
1.) Both gradients were evaluated from temperature dif- 
ferences measured with a pair of Au(Fe0.07%)-Chromel 
differential thermocouples and the distances between the 
thermocouples Ax = 3.55 mm and Ay = 2.0 mm. Since 
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1: (Color online) The magnetic field dependence of (a) 
and (b) the resistivity. The insets of Figs. 1(a) and (b) 

show the magnetic field dependence of -7-^ and the temper- 
ature dependence of the resistivity under zero field, respec- 
tively. 
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the thermal conductivity of chromel is 13.8 W/mK, which 
is comparable to ~ 10 W/mK of bismuth, an error of the 
heat current estimation should be below 0.15% using the 
chromel wire with 50 micron diameter and 2 cm length. 
To avoid leakage of the heat current by radiation and con- 
duction through the air, we attached Cu radiation shield, 
and evacuated the air down to 10 -4 Pa. The loss of the 
heat current was estimated to be below 1%. We observed 
-J- 51 ^ at T due to a small misalignment of the posi- 
tions of the thermocouple. This component was carefully 
subtracted by varying the magnetic field between -3 and 
+3. 

Figure 1(a) shows the magnetic field dependence of the 
temperature gradient -^-- -^ JL decreases with increas- 
ing the magnetic field, which indicates a negative sign 



AT, 
Ay 



at 250 K and 3 T is 

AT, 



of i?TH- The magnitude of 

0.076 K/mm, which is about 29 % of 0.26 K/mm in ^ 
shown in the inset of Fig. 1(a). Thus, -Rth is evaluated 
to be -0.1 T" 1 at 250 K and 3 T. This is consistent 
with i?TH in the reference [lj]. We would like to point 
out that i?TH is -0.39 T" 1 at 0.33 T and 150 K, which 
is one of the largest Rth values among metal elements 



and semiconductors. 



The nonlinear character of 



ATy 

Ay 



is caused by the changes of the physical properties un- 
der magnetic field as shown in Fig. 2. To determine K e \ 
and /in, resistivity p was measured under magnetic field 
as shown in Fig. 1(b). With increasing the magnetic 
field, p remarkably increases. This large positive mag- 
netoresistance is one of the significant properties of the 
bismuth. The inset of Fig. 1(b) shows p(T) under zero 
field. This is also a typical behavior (80 pflcm at 300 K) 
of a bismuth single crystal. 



Using the Fourier's law j x = — n t 



AT X 
ot Ax 



the total ther- 



FIG. 2: (Color online) Temperature dependence of the ther- 
mal conductivity, (b) Hall coefficient, and (c) Hall mobility. 
The inset of Fig. 2(b) shows the magnetic field dependence 
of the Hall voltage. 



mal conductivity Ktot was obtained as a function of tem- 
perature as shown in Fig. 2(a). Ktot is a sum of the lattice 
thermal conductivity K p h and the electron thermal con- 
ductivity K e i (Ktot = K ph + k c \). k c \ is evaluated by the 
Wiedemann- Franz law, k c \ — LqTo (Lq: Sommcrfeld's 
value 2.44 x 10~ 8 Wil/K 2 , a: electrical conductivity). 
At 300 K, the magnitude of n to t is 110 mW/cmK, which 
is close to 100 mW/cmK in the reference [15|]. The lower 
thermal conductivity than that of a conventional metal 
such as Cu, Ag, or Au is attributed to both the relatively 
low electrical conductivity and low K p h- 

Figure 2(b) shows the temperature dependence of the 
Hall coefficient i?n- The sign of i?H is negative, which 
indicates that major carriers are electrons. This is con- 
sistent with the sign of Rth- The magnitude of Rh is 
1 x 10~ 7 m 3 /C at 300 K and 1 T, which shows a low 
carrier concentration n of 6.2 x 10 19 cm -3 . With increas- 
ing magnetic field, i?n slightly increases. This result also 
reproduces Rh measured by Rosenbaum et al. Vc% Using 
the relationship /!jj = Rn/p, the Hall mobility pu was 
derived as shown in figure 2(c). As expected, a rather 
large /in was observed. At around 200 K, peak structure 
was appeared, in which decreased /xh at around room 
temperature is attributed to the increased phonon scat- 
terings of conduction electrons [l?} . 

Now we compare Rth(T) of bismuth with Eqs. (1) 
and (2). The measured k c \, K to t, and /xh in Fig. 2 were 
used for the calculations. Figure 3 shows Rth (T) at 
1 and 3 T. The dotted and broken lines represent the 
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FIG. 3: (Color online) Temperature dependence of thermal 
Hall coefficient at f and 3 T. The dotted lines represent the- 
oretical calculation using a single band model of the thermal 
Hall coefficient. 



calculations. We found that /j,h well reproduces Rth(T), 
which is consistent with the previous results that the con- 
ventional metals obey Eq. (1). On the other hand, Eq. 
(2) that describes Rth(T) of the semiconductor does not 
give a reasonable fitting of Rth(T). The reason why heat 
conduction throuth lattices does not decrease -Rth of bis- 
muth is an open question. Here, we briefly discuss a two- 
band model. As a first order approximation (magnetic 
field B ~ 0), the thermal Hall coefficient -R TH f° r a two 
band model is described as K-R TH — k p Rth p + KnRTHn 
— Kp/iHp — K n fJ.Hn [18], where p and n indexes repre- 
sent a hole- and electron-like band, respectively. In the 
referece [lfj] , a ratio of hole and electron concentrations 
c = — = 0.55, electron mobility /in n = 0.43 m 2 V -1 s _1 , 
and hole mobility fiHp = 0.15 m. 2 V -1 s _1 at 270 K in a 
bismuth thin film were reported. If one applies these pa- 
rameter to the equation shown above, -R TH ^ s evaluated 
to be 93 % of i?THn- This implies that electron trans- 
port is dominant in bismuth. Concerning the deviation 
from Eq. (2), there are two possibile origins; (I) Dirac 
fermion with a linear dispersion, and (2) interband scat- 
tering. Bismuth is known as a Dirac fermion sysytem 
[19J |. where linear-dispersion relation is realized. This re- 
lation gives to* = 0, which can enlarge Rth- Concerning 
the possibility (2), the two band model shown above does 
not take into account the interaction between the elec- 
tron band and the hole band (i.e. interband scattering) 
that explains large diamagnetism in bismuth [20]. This 
effect may increase Rth- As a further analysis, a theoret- 



ical calculation of .Rth on bismuth based on two-bands 
Dirac fermion taking into account interband scatterings 
will be desired. 

Lastly, we would like to add a few comments on an 
improvement of the magnitude of Rh . According to the 
Eqs. (I) and (2), small K p h, large K e \, and large fin are 
essential for large Rth- Mh is described as /zh = m * [„) , 
where e, t, and to* represent electron charge, electron 
scattering time, and effective mass of electron, respec- 
tively. With increasing n, k c \ increases while //h de- 
creases due to increase of to*, which gives an optimum 
value of Rth at an optimum n. To achieve a further 
enhancement of -Rth, one needs to control r. Recent 
theoretical study on thermoelectric performance in bis- 
muth thin film reveals a significant increase of the ther- 
moelectric properties mainly due to the enhanced elec- 
tron scattering time in the Bi thin film with ~ 10-nm 
thickness 21| , in which quantum spin Hall state appears 
due to the strong spin-orbit interaction. To keep a sign 
of the spin momentum in the QSH state, the electron is 
not scattered by impurities, which realizes an increase of 
fin. In the same manner, Rth can largely enhance when 
Hh increases with decreasing a thickness of a bismuth 
thin film. Again, we would like to emphasize that Dirac 
fermion with linear-dispersion relation will give to* = 0, 
which also enlarges Rth- Thus, topological insulators 
and Dirac fermion systems such as Bi2Se3, Bi2Te3 or 
graphene can be a good playground for finding a large 
Rth- The large Rth under low magnetic field in bis- 
muth will open a route for controlling heat current in 
electronic devices. 

In conclusion, we have measured resistivity, Hall coef- 
ficient, thermal Hall coefficient, and thermal conductiv- 
ity of a bismuth single crystal and found that the large 
thermal Hall coefficient is well explained by a theoretical 
calculation. A high Hall mobility and a low thermal con- 
ductivity can cause the large thermal Hall coefficient in 
bismuth. In particular, ~29% of the heat current along 
the bisectrix direction is found to be oriented along the 
binary axis at 250 K and 3 T. We also have discussed that 
the strong spin-orbit interaction can be a driving force of 
a further enhancement of the thermal Hall coefficient. 

We would like to thank S. Murakami for fruitful dis- 
cussion on the thermoelectric properties of Bi thin film, 
and also thank D. Sawaki for technical support on the 
thermal conductivity measurement. 
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